Composite pairing in a mixed valent two channel Anderson model 
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Using a two-channel Anderson model, we develop a theory of composite pairing in the 115 family of 
heavy fermion superconductors that incorporates the effects of /-electron valence fluctuations. Our 
calculations introduce "symplectic Hubbard operators": an extension of the slave boson Hubbard 
operators that preserves both spin rotation and time-reversal symmetry in a large N expansion, 
permitting a unified treatment of anisotropic singlet pairing and valence fluctuations. We find 
that the development of composite pairing in the presence of valence fluctuations manifests itself 
as a phase-coherent mixing of the empty and doubly occupied configurations of the mixed valent 
ion. This effect redistributes the /-electron charge within the unit cell. Our theory predicts a 
sharp superconducting shift in the nuclear quadrupole resonance frequency associated with this 
redistribution. We calculate the magnitude and sign of the predicted shift expected in CeCoIns. 



I. INTRODUCTION 

The 115 family of superconductors continue to attract 
attention for the remarkable rise in the superconducting 
transition temperature from 0.2 K in the parent com- 
pound CeIn 3 [T] to 2.3 K in CeCoIn 5 [2] and finally to 
4.9 K and 18.5 K in the actinide analogues, NpPd 5 Al 2 [3 
and PuCoGa5[3], respectively. This last rise has often 
been attributed to the increasing importance of valence 
fluctuations, and here we seek to make this connection 
explicit. Although the Ce 115s are local moment sys- 
tems, neutron studies of crystal fields indicate broad 
quasi-elastic line-widths comparable to the crystal field 
splitting, indicating valence fluctuations^. The highest 
T c 115s, the actinides, involve 5/ shell electrons, which 
are less localized than their 4/ Ce counterparts, and 
thus are expected to be even more mixed valent. Re- 
cent 237 Np Mossbauer studies of NpAl 2 Pds suggest that 
the valence of the Np ions actually changes as supercon- 
ductivity develops 6 . However, despite different degrees 
of mixed valency, both CeCoIn 5 and NpPdsAl 2 contain 
the same unusual transition from local moment param- 
agnetism directly into the superconducting state, shown 
in Figure [T] 

This direct transition suggests that the localized mo- 
ments play a direct role in the pairing, which previously 
led us to propose t hat the 115 materials are composite 
pair superconductors^^. Composite pair superconduc- 
tivity is a local phenomenon involving the condensation 
of bound states between local moments and conduction 
electrons^ it can be thought of as an intra-atomic ver- 
sion of a d-wave magnetic pairing, between conduction 
electrons in orthogonal screening channels rather than 
on neighboring sites, 



*sc 



CO = Wijv{i<72mj ■ Sfj 



(1) 



Here ip[ 2 j create local Wannier states of conduction elec- 
trons in two orthogonal symmetry channels at site j, and 
Sfj is the local /-moment on the same site. The com- 
posite pair combines a triplet pair of conduction electrons 
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FIG. 1. Local moments are seen in the Curie- Weiss suscep- 
tibilities: CeCoIn 5 (T c = 2.3K ^and NpPd 5 Al 2 (T c = 4.9^ 
are reproduced and rescaled by x{T c ) to show their similarity 
(data below T c not shown). 



with a spin- flip of the local moment, such that the overall 
pair remains a singlet. 

While the close vicinity of the Ce 115s to antiferromag- 
netic order has led to a consensus that pairing in these 
materials is driven by spin fluctuations, the two supercon- 
ducting domes in the Ce(Co,Rh,Ir)In5 phase diagram^ 
suggest that composite pairing may provide a second, 
complementary mechanism in the Ce 115s^. The absence 
of magnetism in the actinide 115 phase diagrams suggests 
that composite pairing may play a more important role. 
Composite pairs are predicted to have unique electro- 
static signature, resulting in a small redistribution of the 
/-electron charge within the unit cell and an associated 
change of the /-valence. Understanding these charge as- 
pects of composite pairing is essential to disentangling 
the relative importance of magnetic and composite pair- 
ing in these compounds. 

These observations motivate us to develop a theory 
of composite pairing incorporating valence fluctuations. 
Thus far, composite pai ring has been studied within a 
two-channel Kondo modeF^, treating only the local spin 
degrees of freedom. In this paper, we study compos- 
ite pairing within a two-channel Anderson model, which 
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permits us to include the charge degrees of freedom and 
model the effects of valence fluctuations on the super- 
conductivity in the 5/ 115 materials. We predict a 
sharp shift in both the /-electron valence and quadrupolc 
moment at the superconducting transition temperature. 
The quadrupole moment shift should manifest as a shift 
either in the nuclear quadrupole resonance (NQR) fre- 
quency, for CeArins or the Mossbauer quadrupole split- 
ting, for NpPdsA^, and we make concrete predictions 
for these shifts in Section |V] 

The two-channel Anderson model is the natural exten- 
sion of the two-channel Kondo model. 

H = H c + J2H a (j)+J2 H ^)- (2) 

j 3 

Here the local atomic Hamiltonian at site j is, 

H a (j) = E Q X Q0 (j) + E 2 X 22 (j) + efX^ij), (3) 

where the X's are the Hubbard operators^ Xqo = 
|0)(0|, X 22 = |2)(2|, and X aa = \cr)(a\. The atomic 
states |0), \a) and |2) are shown in Figure [2j where we 
take the doubly occupied state to be a singlet containing 
/-electrons in two orthogonal channels, 

|2) ee |r 2 ® r x ) s = ^ E *s*(*)fk*fk -*\°)- ( 4 ) 

The orthogonality of these two channels is essential to 
the superconductivity, and we take the Coulomb energy 
for two electrons in the same channel to be infinite, mak- 
ing this model closer to two copies of the infinite-C/ An- 
derson model than to the finite-?/ model 14 . These /- 
electrons hybridize with a bath of conduction electrons, 
H c = J^w a e k c k C r c ko' in two different channels 

Hv(j) =^V 1 4 ja X „(j)+h.c.+V 2 4 ja aX_ a2 (j)+h,c., 

j 

(5) 

where the iprja are Wannier states representing a con- 
duction electron in symmetry T on site j. X$ a — |0)(cr| 
and X 2a = |2)(cr| are the Hubbard operators between 
the singly occupied state and empty or doubly occupied 
states, respectively (see Fig. [3]). Here we have adopted 
the notation of Ce atoms, whose ground state is a 4/ 1 
doublet, but this formalism also applies to Np 5/ 3 atoms, 
where both |0) = |5/ 4 ) and |2) = |5/ 2 ) are J = 4 crys- 
tal field singlets, while \a) represents one of the J = 9/2 
crystal field doublets of 5/ 3 . 

To develop a controlled treatment of superconductiv- 
ity in the two-channel Anderson model, we introduce a 
large- N expansion. Preserving the time-reversal symme- 
tries of spins is essential to capture singlet superconduc- 
tivity in the entire family of large N models, and large 
N methods based on SU(N) spins, like the usual slave 
boson approach, will lose time-reversal symmetry for any 
N > 2. Treating superconductivity in the Kondo model 
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FIG. 2. Virtual charge fluctuations of a two channel Ander- 
son impurity, where the addition and removal of an /-electron 
occur in channels I\ and T2 of different crystal field symme- 
try. The ground state is a Kramer's doublet, while the excited 
states |0), |2) are singlets. The excited doublet \T2 : a) rep- 
resents a higher lying crystal-field level and is excluded from 
the Hilbert space of the problem. The relevant charge states 
of Ce 3+ (4/ 1 ), Np 4+ (5/ 3 ) and Pu 3+ (5/ 5 ) are indicated. 



therefore requires symplectic spins, the generators of the 
SP(N) grourJ 3 , 

^ = /i/ r «wy-«. (6) 

Here N is an even integer and the indices 

a G {±l,±2,...±iV/2} (7) 

are integers running from —N/2 to N/2, excluding zero. 
We employ the notation a = sgn(a). These spin oper- 
ators invert under time-reversal S — > OSd^ 1 = — S and 
generate rotations that commute with the time-reversal 
operator!^. 

Physically, the spin fluctuations of a local moment are 
generated by valence fluctuations. Theoretically, these 
valence fluctuations can be described by Hubbard opera- 
tors, and in a symplectic- generalization of the Ander- 
son model, anticommuting two such Hubbard operators 
must generate a symplectic spin, satisfying the relations: 

{X 0a , Xp } = X a p + X m 5 a p (8) 
= S a p + ^X 00 H — ~j^J $a,p, 

where the last equality follows from the traceless defi- 
nition of the symplectic spin operator, S a p — X a p — 
^l^-Sap. In Section |hJ we show that a proper symplec- 
tic representation requires the introduction of two slave 
bosons to treat the Hubbard operators in a single chan- 
nel: 

X 0a = f a + a*af]_ a 

X 00 = tfb + aV (9) 
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As we shall demonstrate, these symplectic Hubbard op- 
erators maintain the invariance of the Anderson Hamilto- 
nian with respect to the SU (2) particle-hole transforma- 
tion, f a — > cos 6* f a + sign(cr) sin^/lj^. This extension 
of the slave boson representation with local SU{2) gauge 
symmetry was originally introduced by Wen and Lee as a 
mean-field treatment of the t — J modef^. Here we show 
that these Hubbard operators maintain the SU{2) gauge 
symmetry for arbitrary SP(N) versions of the infinite U 
Anderson and t — J models. This SU(2) gauge symmetry 
is essential for eliminating the false appearance of s-wave 
superconductivity in the finite-t/ Anderson model, where 
the two channels are identical. 

What do we learn from this Anderson model picture 
of composite pairing? A key result is that the amount 
of composite pairing can be written in a simple gauge- 
invariant form as 



Asc = (A 02 ) oc (ip\a (ia 2 ) V4 



Sfh 



(10) 



where the Hubbard operator Xq2 — |0)(2| mixes the 
empty and doubly occupied states. This mixing resem- 
bles an intra-atomic version of the negative-?/ pairing 
state, but double occupancy costs a large, positive en- 
ergy. The physical consequence of this mixing is a redis- 
tribution of the charge in the /-electron orbitals, as the 
three charge states: empty, singly and doubly occupied, 
all have different charge distributions. Such a rearrange- 
ment not only changes the total charge in the /-shell, 
a monopole effect, but also will result in a quadrupolc 
moment associated with the superconductivity. This 
quadrupolc moment will lead to modified electric field 
gradients that should be detectable as a sharp shift in the 
nuclear quadrupolar resonance (NQR) frequency at the 
neighboring nuclei, or as a sharp shift in the Mossbauer 
quadrupolar splitting at the /-electron nuclei. 

This paper is organized as follows. First we intro- 
duce the symplectic Hubbard operators in section [H] and 
demonstrate that the symplectic- A and SU(N) large- 
N limits are identical for a single channel. We then 
generalize this formalism to the case of two channels 



in Section III and show how composite pairing natu- 



rally appears as a mixing of the empty and doubly oc- 
cupied states. The mean-field solution is presented in 
Section |TV] in which we demonstrate how the supercon- 
ducting transition temperature increases with increasing 
mixed valence. In Section [V] we calculate the charge dis- 
tribution of the /-orbitals in the state with composite 
pairing and make a concrete prediction for a shift in the 
NQR frequency at T c in CeCoIns . Finally, section VI dis- 



cusses the implications for the finite-C/ Anderson model 
and examines the broader implications of our results. 



II. SYMPLECTIC HUBBARD OPERATORS 

Composite pairing was originally discussed in the 
two channel Kondo modeP^, where it is found in the 
symplectic- iV limiP, which maintains the time-inversion 



properties of spins in the large N limit by using sym- 
plectic spins. Correctly including time-reversal symme- 
try allows the formation of Cooper pairs, and thus super- 
conductivity. In order to treat composite pairing within 
a large N Anderson model treatment, we would like to 
develop a set of Hubbard operators that maintain this 
time-reversal property in the large N limit. Hubbard 
operators, like X a {, = \a)(b\, are projectors between the 
local states, |0,ct) of the infinite U Anderson model, de- 
scribing both charge, Xo a and spin, X acr > fluctuations. 

Starting from a spin state \a), hopping an electron off 
and back onto the site generates a spin flip: this condition 
defines the Hubbard operators within a single channel, 
where the Hubbard operators must satisfy a graded Lie 
algebra in which the projected hopping operators anti- 
commute, satisfying the algebra 

{X a0 , Xop} = X a p + X 00 S a j3 (11) 

Now the traceless part of X a p defines a spin operator, 

G — Y ^<*P V 
Sufi — -*~af) — -A 77 

so that quite generally, Hubbard operators must satisfy 
an algebra of the form 



{A q q, Xqp) — S a p 



X (l 



X 
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lap 



(12) 



Thus the commutation algebra of the Hubbard operators 
generates the spin operators of the local moments. In the 
traditional slave boson approach^, the Hubbard opera- 
tors are written X$ a = W f a . The spin operators gen- 
erated by this procedure are the generators of SU(N), 
since 



{fib, tffp) = (flfp - + (fttft + ^) S a0 (13) 

where is the well-known form of SU (N) spins. It 

is thus no wonder that this approach cannot treat super- 
conductivity in the large N limit, for there are no spin 
singlet Cooper pairs for SU(N) with N > 2. We require 
instead that the spin fluctuations generated by Hubbard 
operators are symplectic spin operators 



(14) 



We now show that the Hubbard algebraQ can be sat- 
isfied with symplectic spins, using the introduction of two 
slave bosons, 



X a0 = fib + &f- a a, 
X Q p = b^f l3 + o^/L 



(15) 



These operators satisfy the Hubbard operator algebra 
(111, but now 



X, 



a/3 



fifp + &Pf- a fL - s aP + s, 



a/3 



(16) 
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describes a symplectic spin operator while 
X 00 = tfb + a f a 



(17) 



is the representation of the empty state operator. 

At first sight, the expedience of this new representa- 
tion might be questioned: why exchange the simplicity of 
the original Hubbard operators, X a b — \a)(b\ for a pro- 
fusion of slave boson fields? However, while the original 
Hubbard operators may appear to be simple, they are sin- 
gularly awkward to treat in many-body theory^, due to 
their noncanonical anti-commutation algebra, (p]). The 
slave boson representation allows us to represent these 
operators in terms of canonical bosons and fermions. 
By doubling the number of slave bosons the symplec- 
tic character of the spins is preserved at all even val- 
ues of N, and we shall see that this process encodes the 
hard-to-enforce Gutzwiller projection as a mathemati- 
cally tractable SU (2) gauge symmetry. 

Physical spins are neutral, and thus possess a continu- 
ous particle-hole symmetry. This property is maintained 
by symplectic spin operators, which can be seen most 
naturally by introducing the generalized pair creation op- 
erators, 



(18) 



Q>0 



which allow us to construct the isospin vector, ^ 
* 2 , ^3), where 



#) , * 2 = -i - 
*3 = J2(flfa - f-afL) =n s - N/2, (19) 



Q>0 



and n s = ^2 a (&\f a f a \ a ) is the fermion number in the 
singly-occupied states, i.e. n s plays the role of the num- 
ber of spins in the above constraint vf^ . The isospin vec- 



0, show- 



tor commutes with symplectic spins, 

ing that the symplectic spins possess an 5/7(2) gauge 
symmetry: a continuous particle-hole symmetry that al- 
lows us to redefine the spinon, f a — > uf a + vaf^_ a . This 
symmetry is reflected in the requirement of two types of 
bosons, as the empty state does not distinguish between 
zero and two fermions, and thus requires two bosons to 
keep track of the two ways of representing the empty 
state, &t|0) and at/|/J|f2) = a^^\Q), where |fi) is 
the slave-boson vacuum. Of course, there is only one 
physical empty state, as becomes clear when we restrict 
these Hubbard operators to the physical subspace. In 
order to faithfully represent the symplectic spins, the 
sum of the spin and charge fluctuations must be fixed, 
S 2 + # 2 = f (y + 2)0. While this constraint is en- 
forced by setting = in the pure spin model, here 
we must equate our two types of charge fluctuations, by 
setting 



a) a 



Q>0 

Q- = J2f-^f a + ^b = o. 



(20) 



Q>0 



The three operators (Q±,Qs) commute with the 
Hamiltonian, which imposes the constraints associated 
with the physical Hilbert space. The constraint reflects 
the neutrality of the spins under charge conjugation: 
Q3 conserves total electromagnetic charge, and prevents 
double occupancy, while Q± kills any states with s-wave 
pairs on-site. From the form of Q3, we see that b and a 
have opposite gauge charges. The gauge invariant states 
satisfying the constraint (for N = 2) are, 



|0) = (fct + a t^t) |fi). 



(21) 



The symplectic Hubbard operators can be written 
more compactly by using Nambu notation, 



where 



XQ a — f a , Xq — B^B 



(22) 



(23) 



The constraint becomes Q = B^fB + J2 a >o fa? fa = 0- 
For N = 2, these Hubbard operators are the SU(2) 
slave bosons introduced by Wen and Lee in the context 
of the t — J model 16 . Here, it becomes clear that the 
SU(2) structure is a consequence of symplectic symme- 
try, present in both the symplectic- N Kondo and Ander- 
son models, not just for N — 2, but for all N. This sym- 
metry can be physically interpreted as the result of va- 
lence fluctuations in the presence of a particle-hole sym- 
metric spin. 

The Hamiltonian of the one channel, single impurity 
infinite-?/ Anderson model can now be expressed via sym- 
plectic slave bosons as follows: 

H = J2 £ kcL c ka + E {tfb + a} a) 

k . a 

+ J2vi>i(b^f a +a^fl a ) + h.c. (24) 

a 

The SU(2) gauge symmetry of this model becomes 
particularly evident in the field-theoretical formulation, 
where the partition function is written as a path integral: 



Z= / V[ft,$,pJ,Bi,B]t 







(25) 



where the Euclidean action S is written in Nambu nota- 
tion as follows 
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( 

S = f dT £ ^L(^+e fc r 3 )4a+^ i ^Tr [r 3 ^t(a T + E + \j ■ t)A,] + ]T [f] a (d T + X j ■ f)f ja + (v$ a A]f ja + h.c 



(27) 



where ip^ — (ipa> Qi^-a)) f3 = 1/T is the inverse temper- 
ature and the SU (2) matrix Aj of slave bosons is defined 
as 



i 



(28) 



The first three terms in Eq. (271 describe the dynamics 



of conduction electrons ip, /-electrons / and the slave 
bosons B respectively, with the constraint Q = imposed 
by introducing the vector of Lagrange multipliers A. The 
last term describes the coupling between the f-atom and 
conduction electrons. 

The action (27) is manifestly invariant under the 



SU(2) gauge transformation via unitary matrix gf 



Z ■ f 



► A,- • r 



(29) 



At first sight, the a hybridization term in the Hamilto- 
nian, (24) appears to give rise to (c'f') pairing terms, 



however, the SU(2) symmetry allows to redefine the 
fields, making the substitution 



B = 



gB = 



(30) 



This transformation eliminates the composite s-wave 
pairing, f a + at f_ a a — > b^f' a , recovering the usual 
U(l) slave boson Hamiltonian. 

One of the the most important physical consequences 
encoded in this local gauge symmetry is the complete 
suppression of any s-wave pairing in the single channel 
Anderson model. Indeed, the constraint operator Q + = 
S Q >o faf-a ma y De directly interpreted as the on- 
site s-wave pair creation operator plus a term wa that 
can be removed by the above gauge transformation. In 
particular, the mean-field constraint 



(Oh 



(s-wave pair creation operator) = 



hard-wires the strong suppression of local s-wave pair- 
ing into the formalism. In the language of superconduc- 
tivity, the constraint Q± — plays the role of an infi- 
nite Coulomb pseudopotential, or renormalized electron- 
electron interaction, /i* . It is the satisfaction of this con- 
straint that drives anisotropic composite pairing. 



III. THE TWO CHANNEL ANDERSON 
LATTICE MODEL 

We turn now to the two channel lattice Anderson 
model, where the two channels involve charge fluctua- 
tions to the "empty" and "doubly" occupied states; here 



we use the language of Ce 3+ whose ground state config- 
uration is / , but this model captures any system with 
valence fluctuations from f n f n , where n is 

odd. For n = 1, the empty state, |0) is trivially a singlet, 
and we choose the doubly occupied state to be a singlet 
formed from electrons in two orthogonal channels, 



|2> = 5/k >a &,_ Q |0>, 



(31) 



where |Fi : ±) = |±) is the ground state crystal field 
doublet, and |T2 : ±) is an excited crystal field doublet. 
In the finite-?/ Anderson model, the symmetries of the 
electron addition and removal processes are identical, but 
here Hund's rules force the second electron to be placed 
in a channel orthogonal to the first. 

The local atomic Hamiltonian at site j can be ex- 
pressed in terms of Hubbard operators, 

H a (j) = E X 00 (j) + E 2 X 22 (j) + J e/ Z„0) 

+ ^ A CT l r2(; . r2ff , (32) 

IT 

where the X's are the Hubbard operators projecting 
into the empty, Aqo = |0)(0| and doubly occupied, 
A22 = |2)(2| states, as well as the projectors into the 
ground state, X aa = |r 1 a-)(r 1 cr| = |ct)(ct|, and excited, 
Xr 2 <7:r 2 <T — |r20-)(r2<r| crystal field doublets. If we mea- 
sure the energies from the /-electron level, Eq — — e/ and 
E 2 = U\ 2 + £/ are both positive, where Ui 2 is the Hub- 
bard U for the doubly occupied state with one electron in 
Ti and the other in T 2 . We take U\\, U 22 — > 00. The /- 
electrons hybridize with a bath of conduction electrons, 
H c = X)k a £kC kcr Ckcr m t w0 different channels, 

H v (j)=J2Vi4 j *XoA3)+V2*4 j(T X- (r2 (j)+h. C . (33) 



The angular momentum dependence is hidden inside the 
conduction Wannier states, 



ifcjc 



rkj 



Cko 



(r = l,2) (34) 



where the crystal field form-factors [3>rkL CT ' = {kY(j'\ka) 
are proportional to unitary matrices. X( )a = \0)(Tia\ 
and X 2a = |2)(ricr| are the Hubbard operators between 
the singly occupied ground state and empty or doubly 
occupied states, respectively. 

The symplectic--ZV Hubbard operators describing 
charge fluctuations in the two channels are a simple gen- 
eralization of those used for a single channel. For a single 
site, these are given by 



Xoa 



b\f a + %afL a , X Q = &0&0 + <4 a 0; 



G 



X 2a = 4/a ~ 02«/-a) ^22 = 4 6 2 + 4 a 2: (35) 

which can again be written more compactly by using a 
Nambu notation, 



Xo a — -Bq/q, X 2a — B 2 f a 

Xqq — BqBq, X22 = B\B 2 , 



where 



bo 
«o 



Bo 



b 2 

-a 2 



(36) 



(37) 



We are free to choose the sign of a 2 , and in the above 
we have chosen the negative sign to preserve continuity 
with the results from the two channel Kondo model 7 . 
The doubly occupied state is then, 



\2) = (bt-4&)\n) 



(38) 



where ^ = J2 a >o faf-a i s the pair creation operator, 
as before. 

In the two channel model case the constraint becomes, 
Q = BlfB Q + B\rB 2 + P a Tf a = 0. (39) 

The intersection of the two Hubbard algebras gives rise 
to an extra doublet, 



bl4 + bWo)fl\n), 



(40) 



which we interpret as the excited crystal field doublet 
because it is reached by destroying a Ti electron from 
the doubly occupied state, leaving behind the T 2 elec- 
tron. However, this extra doublet is not killed by X 00 
or X 2 2, which are no longer the projectors within this 
larger space. The full Hamiltonian term for this state 
must therefore be quartic in the bosons, and subtract off 
the contribution from X 00 , X 22 and X aa , 



the ordered state proposed for URu9Si9|19|. which mixes 
ground state and excited singlets. 

The operator X 2 o = |2) (0| mixes the empty and doubly 
occupied states, and can be thought of as a pair creation 
operator. As we show in the next section, this opera- 
tor acquires an expectation value when composite pair- 
ing develops. The development of an intra-atomic order 
parameter (X02) is reminiscent of pairing in a negative 
(attractive)-?/ atonPHl, but unlike a negative-17 atom, the 
empty and doubly occupied sites are excited states of the 
atomic Hamiltonian, and they are only partially occu- 
pied as a result of valence fluctuations. Furthermore, the 
paired state is an anisotropic singlet with nodes. In the 
case of the Ce 115 materials, this product has d- wave 
symmetry^. 

Again, we should consider the question of supressing 
s-wave superconductivity. Here, 

Q+ = E tif-« + b o«o + bW 

We know from the one channel model that ao can be re- 
moved by a gauge transformation, eliminating the mid- 
dle term, but the b\a,2 term cannot be uniformly elim- 
inated. However, it will only be nonzero in a state 
where both (X 2 o) and (^ria;r 2 /3) are nonzero, imply- 
ing coexisting superconducting and composite nematic 
order. Under normal circumstances, these two phases re- 
pel one another, and s-wave superconductivity is wholly 
suppressed. 



IV. THE LARGE N TWO CHANNEL 
ANDERSON LATTICE MODEL 

We are now able to write the symplectic-./V two chan- 
nel Anderson lattice model in terms of the symplectic 
Hubbard operators, 



{A CF -(E + E 2 +e f )} (bl4 + albl) (b a 2 + a b 2 ) , (41) H = £ + £ Eq bI(j)B q (j) + E 2 B\{j)B 2 {j) 



where Aqf is the crystal field effect splitting. We neglect 
this term for much of our analysis, as it turns out to be 
irrelevant for the superconducting states of interest. 

In addition to an extra state, there are two additional 
operators, 



Vi 



Vo 



3 

E4^/!-a5 2 (j)+h.C. 



(43) 



A20 — {X2a, X a o} — BIBq — b 2 bo — a\ao 



X 



r lQ; r 2 /3 



= -{X 0a ,X 2 p}=Bl(-iT 2 ) 



Bl 



(42) 



The operator ^ria;r 2 ^ mixes the ground state and ex- 
cited crystal field doublets, leading to a composite den- 
sity wave state, (cj//^). As this state mixes two ir- 
reducible representations of the point group, it necessar- 
ily breaks the crystal symmetry and this phase can also 
be called a composite nematic. This mixing resembles 



In order to keep the Hamiltonian extensive in N, we have 
rescaled the hybridization terms by (7V/2) -1 / 2 , so that 
they recover the correct N = 2 form, ( |33| . This rescaling 
implicitly assumes that the slave bosons fields Bi (I — 
0,2) will acquire a magntitude of order O(VN). Since 
there are only two flavors of bosons, for the bosons to 
play a role in the large N limit, they must condense. 

Following the construction that led to the path integral 
description of the single-channel Anderson model (27) 



above, one can write down the Euclidean action that cor- 
responds to this Hamiltonian: 
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FIG. 3. The six physical states in the Hilbert space of the N — 2 two channel Anderson model, representing the ground state 
and excited crystal field doublets, and the empty and doubly occupied singlets. Arrows show the operators that move between 
states in the Hilbert space. 



P 

S 2 [i>J,A ,A 2 ] = dr $kJP* + £ fc T 3)V'ka 

k.a>0 
3 v 



TsAl^dr + E + X r r)A 0j 
Vi 



-Tr 



2.; 



(44) 



where ^\ja = Wrja' s S a ( a )' l Prj-a) are Nambu spinors 
for the Wannier states in each channel. Above, we have 
collected the slave bosons into the SU(2) matrices, 



A ] - 



b\ aj 
a -b 



v4q 



(to 



(45) 



The action (44) remains manifestly invariant under the 
SU(2) gauge symmetry B\ — > cjjBi, where c/j is an SU{2) 
matrix. The matrices Ao, A2 transform under this gauge 
symmetry as At — > gjAi, leaving the product 



b a 2 + b 2 a a\a 2 - b J b 2 

t„_ j,t„t 




b% 



a 2 a,o b 2 a' n 



atb} 



(46) 



gauge invariant. 

The off-diagonal components of the product .A^O 
have the physical meaning of composite pairing between 
conduction and /-electrons. Indeed, in the Kondo limit, 
Eq,E 2 3> npV 2 , we can connect the two channel An- 
derson model results to those from the two channel 
Kondo model 7 . A Schrieffer- Wolff transformation maps 
bo — > Vi, ap — > Ai, b 2 — > A 2 and a 2 — > V 2 , where Vp and 



Ar are the particle-hole and particle-particle hybridiza- 
tions, respectively. We can now identify the off-diagonal 
components of A\Aq oc ViA 2 — V 2 Ai with composite 
pairing. More generally, the components of can be 

identified with the state mixing operators, 



A\A a 



-^20 



— X02 
X T 2 T 1 



(47) 



which confirms the identification of (X02) with composite 
pairing, and implies that the composite pair state con- 
tains an admixture of the empty and doubly occupied 
states. 

To write down a translationally invariant Hamiltonian, 
we assume that the expectation values of the slave bosons 
are uniform, which allows us to write down the Hamil- 
tonian in momentum space. To simplify this step, we 
temporarily drop the spin-orbit dependence of the Wan- 
nier functions, treating the form factors as spin-diagonal, 
Tprja = Cka^rke -11 ' ^ ■ This allows us to absorb the 
momentum dependence of the form-factors into the hy- 
bridizations by defining Vrk = - r=^<t>vk- To obtain the 



d-wave symmetry which arises naturally from the spin- 



orbit form factors^, now we must explicitly make V\kV 2 k 
d-wave. The full spin-orbit dependence can be restored 
in a similar manner to the two channel Kondo model 7 . 
The Hamiltonian is, 



H =J2(^cl a c ka + V lk cl a Blf ka + V 2k cl a aft_ a B 2 + h 

ka 

hK (e bIb + E 2 B\B 2 ) , (48) 



where Af s is the number of lattice sites. As before, the 
SU(2) constraint is implemented with a vector of La- 
grange multipliers, A, 



A ■ 



B^tBo 



B 2 tB 2 



(49) 



where we have also assumed that A is translationally in- 
variant, enforcing the constraint on average. This ap- 
proximation becomes exact in the large N limit. 

The Hamiltonian can be re-written in the compact 
form with Nambu spinors, 



#=E( 5 L/L 



ka 



e k r 3 V Lk .J&+V 2k J&fc k , 
Vi k A +V 2k A 2 A-f 



[E BlB + E 2 B\B 2 + A • (b\tB + b\tB, 



,(50) 



We note that while in the one-channel infinite-?/ An- 



derson model (27), the spurious superconducting s-wave 



state was eliminated by fixing the SU(2) gauge transfor- 
mation, the two-channel model (44) generally possesses 



a composite superconducting ground state, which can- 
not be eliminated by a gauge fixing procedure. However 
we can always define a transformation g such that do 



vanishes, so that B^ — > (6 T ,0). In this new basis, the 



composite pairing operator Xq 2 = B\B 2 (42) takes on 



the simple form, Xq 



A. The mean field solution 

In order to study the effects of mixed valence on the 
superconductivity, we examine the mean field solution of 
the composite pair state in the symplectic-A^ limit. The 
bosons are replaced by their expectation values, Bi — > 
(Bi) ~ 0(y/(N)), and the Hamiltonian |50| ) becomes 
quadratic in the fermions, which may be integrated out 
exactly, leading to the effective action S e f / [B , B 2 ] in 
terms of the boson fields only. The mean-field solution is 
then obtained as a saddle-point of this action: 



5S e ff[B , B 2 ] 
SB, 



0. 



(51) 



We use the SU (2) gauge symmetry to eliminate the ao 
boson, and now the composite pair state is defined by the 
nonzero expectation value of {X a2 ) = {b\b 2 ). While the 
a 2 boson can in principle acquire an expectation value, 



it would lead to a uniform composite density wave so- 
lution, which is generally unstable to the composite pair 
solution. Therefore in what follows, we shall set (02) = 0. 
The resulting free energy can be rewritten in terms of our 
mean field parameters, where we replace (6;) / 'y/~N/2 with 
bi for clarity, (however keep in mind that these have lost 
their dynamics), 

F = -NT £ log 2 cosh 

k± 

+ I ^l[bl(E + X 3 )+bl(E 2 + X 3 )]. (52) 



J\f s is the number of sites, and the dispersion of the heavy 
electrons is given by four branches: U3).± and — Wfc± , where 
uj k± = \/a k ± r fc , and 



a k =b 2 + + -(el + \l + \l), T k 



ll = [e k A 3 - b 2 _] + [2V lk b V 2k b 2 - ekAi] 2 . (53) 
We have also defined 



In a nodal composite pair superconductor, the A con- 
straint reduces to A3, as the Ai constraint acts as a 
Coulomb pseudo-potentiaP^ eliminating s-wave pairing, 
and it is thus unnecessary when we choose VikV2k such 
as to give nodal superconductivity. However, if we were 
to treat the finite U model, where V\k = V 2 k, this con- 
straint is essential to eliminate the appearance of a false 
s-wave superconducting phase. 

The mean field parameters are determined by minimiz- 
ing the free energy with respect to 60, b 2l A3 and Ai- To 
understand their implications, we first present the mean 
field equations in real space, 



E 
E 2 



(bij 

a 



N/2-(b Q3 ) 2 -(b 2j ) 2 



(54) 



These equations bear a strong resemblance to the two 
channel Kondo equations, where (60) plays the role of 
the hybridization V± in channel one, while (62) plays 
the role of the pairing field A2 in channel two; here 
the hybridizations are explicitly identified as the magni- 
tude of the valence fluctuations to the empty and doubly 
occupied states. The number of singly-occupied levels, 
n s = J2 a fafa, is no longer fixed to A^/2 and instead 
decreases as hybridization and pairing develop. 

To calculate the phase diagram, we return to the mo- 
mentum space picture, and derive the three equations rel- 
evant for composite pair superconductivity with a nodal 
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order parameter. The last two equations in (54 1 impose 
the constraint Q = 0, fixing n s = N/2 — 6§ — fc»2 an d an- 
nihilating any s-wave pairs, while the first two equations 
determine the magnitude of the valence fluctuations to 
the empty and doubly occupied states, respectively: 




where 



.4 



X 3 a k + e k \b 2 _ - A 3 e k ] \ 
Aia k + e k [2F lk F 2k &o&2 - ekAi] 



(e k + A 3 ) 2 + Af 
(e k - A3) 2 + Xi 



2Ail4kV A 2 k e k 6 2 
2AiVikV2 k e k 6o/ 



.(56) 



In addition to these four equations, we must also fix the 
total electromagnetic charge in the system by keeping 
the total number of conduction electrons plus physical 
/-electrons constant. Notice that the number of physical 
/-electrons, nf = N/2 — b\ J vb\ counts the electrons in 
both the singly and doubly occupied states, and differs 
from the occupation of the spin states, n s = N/2 — b\ — b\ 
by 2b\. 

These equations can be solved numerically for a sim- 
ple two-dimensional model where we take an s-wave 
Vik = V and a d-wave V 2 k = V (cos k x — cos k y ), and 
take the two-dimensional conduction electron dispersion, 
£k = -2i(cosfc 

x ~\~ cos ky) /i, where fi is adjusted to fix 
the total charge, nt +n c . Ai can be set to zero here since 
by construction, we are not considering s-wave composite 
pairing. This model contains three non-trivial phases: 



0.15 



0.05 




n 2 0.03 



FIG. 4. The superconducting transition temperature is 
plotted against the occupancy of the empty, no and dou- 
bly occupied, 712 states, showing that T c typically increases 
with increasing mixed valency. The black line indicates the 
maximal T c , where T* = T£ , and its curve is due to the 
d-wave nature of the second channel. Two different phys- 
ical paths are shown: (1) The cost of doubly occupancy, 
U12 = E% — Eg — 2 1 e/ 1 is tuned while fixing the /-level, e/ 
(blue, solid path). Here, Eo is fixed and E2 increases from 
left to right. (2) Alternatively, the hybridization, V can be 
increased, as by pressure or by exchanging 4/ for 5/ ions 
(green, dashed path). T c is scaled by A = tyVi, and we have 
fixed n c = 1.8 in units of N/2. 



always finite and the ground state will always be 
superconducting. The /-electron valence will gen- 
erally differ from N/2, with nonzero contributions 
of both /° and f 2 . 



• For Eq -C E 2 , a heavy Fermi liquid develops be- 
low the temperature T* . In this phase, 60 becomes 
nonzero and the hybridization has symmetry T\ 
due to screening of /-spins by conduction electrons 
in this channel. The /-electron valence will be less 
than N/2, and these results will be identical to the 
infinite- U Anderson model discussed in the intro- 
duction. 



For E 2 < E 



T 2 , where 62 becomes nonzero 



a heavy Fermi liquid develops at 
This Fermi liquid 
will have the symmetry of the excited doublet, T 2 
and the /-electron valence will be larger than N/2. 
Again, these results should be identical to the ap- 
propriate infinite-?/ Anderson model. 

• Below the superconducting temperature, T c , a 
phase with d-wave composite pairing develops out 
of the heavy Fermi liquid, where becomes 
nonzero. T c < {Tf,T£} is maximal where T c = 
T* = T 2 * and here the composite pair superconduc- 
tor develops directly out of the high-temperature 
state of free spins, bypassing the heavy Fermi liq- 
uid. As superconductivity is driven by the Cooper 
channel in the heavy electron normal states, T c is 



Note that in the full model, the Kondo temperatures Tj* 
and T 2 * mark crossovers into the heavy Fermi liquids. The 
appearance of phase transitions associated with conden- 
sation of (61) and (62) here is a spurious consequence of 
the mean-field large-A^ treatment and is resolved with 
1/iV corrections^. However, the superconducting phase 
transition is not spurious and will survive to finite TV. 

In Figure [4j we plot the superconducting transition 
temperature versus the proportion of empty states, no = 
Oq and that of doubly occupied states, n 2 — b\. Larger 
uq and n 2 indicate a greater degree of mixed valency, 
which may be obtained by varying V, Eq and E 2 . We 
show two possible paths for tuning real materials: (1) 
By varying U\ 2 , E 2 = et + U12 will change, while Eq re- 
mains fixed. This path is qualitatively identical to that 
of the two channel Kondo model when J2/J1 is varied: 
there is a superconducting dome with maximal T c where 
T* = T 2 *; (2) By changing the hybridization, V and keep- 
ing Eq/E 2 fixed, T c can increase monotonically with in- 
creasingly mixed valence, as it does between CeCoIn5 
and PuCoGas. A similar increase slightly away from the 
maximal T c can explain the non-monotonic change in T c 
with increasing pressure in CeCoIn5 31 . 
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V. CHARGE REDISTRIBUTION 

As the development of composite pairing mixes the 
empty and doubly occupied states, each develops a non- 
zero occupation, and the charge density changes. The 
link between the /-electron charge and the development 
of the Kondo effect in the one-channel Anderson model 
was explored by Gunnarsson and Schoenhammer^l, who 
showed that the /-electron valence, nj decreases grad- 
ually with temperature through the Kondo crossover, 
n f (T) = 1 — 6 2 (T) ■ The mixing of the empty and doubly 
occupied states adds a new element to this relationship, 
and the consideration of real, non-s-wave hybridizations 
allows us to explore the higher angular momentum com- 
ponents of the charge distribution. The charge density 
can be written, p(x) = -0j l (x)'0 Q ,(x), where f/> a (x) creates 
a physical electron of spin a at x. The electron field x 
can be approximately decomposed as a superposition of 
two orbitals Ti and T2 at nearby lattice sites j, 

(57) 

where we have reintroduced the spin-orbit form factor, 
[^r]a/3 in order to model real materials. The charge den- 
sity of an /-electron located at the origin in channel T is 
p r (x) = Tr|<I> r (x)| 2 i?(x) 2 , where i?(x) is the radial func- 
tion for the /-electron, and |$r(x)| 2 is a diagonal matrix. 
If we assume that the overlap of electrons at neighbor- 
ing sites is negligible, the total charge density has three 
different terms, 



<!>■; 



a/3 



(x-R.O/L^+h.c, 



a/3 



(58) 



as it 



where we have kept the spin indices of 

may not be diagonal in spin space. In terms of the Hub- 
bard operators, we can replace, 



fhpfw = X M + n *(j) = BljB 2j + n s (j) 
fljpfw = *M + Xr a pr a ()U) = b\ } B 23 + 7 t 7 
fl ja h j0 = X riar20 (j), (59) 

where as before, the number of singly occupied states 
on site j, n s (j) = N/2 - bLBqj - BLBy, and 7 t = 



-Bq(— i<72)[B\] T — 6g a 2 + ^2 a o * s ^ ne operator from equa- 



«.t„t 



ht„t 



tion (41|, so that 7'/! creates an electron in the ex- 
cited single-electron crystal field state ll^a), see Fig. [3] 
The third term, Xr lQ; r 2( 3 = a)(r2, /3| mixes the two 
crystal field states, shifting charge from Ti to T2. In 
the large- N mean-field theory, this state corresponds to 
non- vanishing diagonal elements of the matrix (A\Ao) in 
Eq. (46), (A ri ar 2/ 3) ~ ao&2 + 02^0- Since Ti and T 2 are 



two different representations of the crystal point group, 
such mixing implies that the crystal symmetry is sponta- 
neously broken: in Ce-115 materials it would correspond 



to orthorhombic "nematic" density wave state, which has 
not been observed. We are primarily interested in super- 
conducting instability, and therefore ignore such a state 
by setting a 2 = in what follows, assuming we have al- 
ready gauge fixed clq — 0. Any occupation of the excited 
crystal field state, (^r 2 ^r 2 /3) is a l so eliminated by this 
Ansatz. 

In the heavy Fermi liquid and superconducting states, 
we may use the constraint to rewrite fljpfijp = N/2 — 
Nao(j). The charge distribution becomes, 



p(x) =^pi(x-Rj) 



N 



-Xoo(j) 



^2(x-R,)A2 2 (jX-60) 



Integrating this charge density around a single site gives 
us the /-electron valence, nt = 1 — 6§ + 6| , which we plot 
as a function of T in Figure [5] for two different E2/E0. 
The phase transition at T* is an artifact of the large iV 
limit and becomes a crossover for any finite N, but the 
sharp kink in nf at the superconducting transition tem- 
perature remains for all N. Experimentally determining 
the /-electron valence may be possible using core-level 
valence spectroscopy for Ce compounds or the Mossbauer 
isomer shift for NpPdsA^, where current measurements 
do indicate a small, positive change in the isomer shift 
through T c 6 . If Np is in the 5/ 3 valence state, with 
dominant fluctuations between 5/ 3 ^ 5/ 4 , then the iso- 
mer shift will increase with decreasing temperature down 
to T c , where it will begin to decrease sharply as 5/ 2 states 
mix in, as the black curve in Figure[5]shows. Experimen- 
tally, the isomer shift was measured at 10K, well above 
T c and then below T c , so the observed positive shift could 
just be due to the increase above T c , and further mea- 
surements are necessary. The clear observation of a sharp 
negative shift precisely at the superconducting transition 
temperature would indicate the presence of composite 
pairing. 

As superconductivity develops, the occupation of the 
doubly occupied state acquires an expectation value, 
leading to an increase in the T 2 charge density. This 
redistribution of the /-electron charge results in a 
quadrupole moment associated with superconductivity. 
Again, since the development of superconductivity is 
a phase transition, the quadrupole moment changes 
sharply at T c . The quadrupole charge component has 
an indirect effect on the superconducting transition tem- 
perature through its linear coupling to strain, leading to 
a linear dependence of T c on the tetragonal strain, c/a. 
Such a linear increase of T c with c/a has been observed in 
both the Ce and Pu 115s, although it is conventionally 
attributed to dimensionality effects^. The quadrupole 
moment of the composite condensate provides an alter- 
nate explanation, and, in addition, leads to the develop- 
ment of electric field gradients around the /-lattice sites, 
which can also be measured directly as a shift in the 
nuclear quadrupole resonance frequency, Ai/jvqk at the 
nuclei of the nearby atoms. 

To make contact with potential experiments, we exam- 
ine CeCoIns in more detail. 115 In atoms have a nuclear 
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CeCoIn 5 


In(l) V zz 


In(2) V„ 


1 7 


-1.5g 7+ ■ 10 20 V/m 2 


+5.0g 7+ • 10 19 V/m 2 


r 6 


+lAq 6 ■ 10 20 V/m 2 


-3.6g 6 • 10 19 V/m 2 



TABLE I. Estimated electronic field gradients at the two In 
sites in CeCoIns due to egr charge in the /-electron orbital, 
r = Iy , Tq, where e is the charge of an electron. 



The NQR frequency measures the electric field gradi- 
ents at two different indium sites in the crystal: the in- 
plane, high symmetry In(l) sites, which sit in the center 
of a square of Ce atoms, and the lower symmetry out-of- 
plane In(2) sites, which are above and in-between two Ce 
atoms. The NQR frequency is given by^, 



FIG. 5. The /-electron valence, nj changes as the Kondo 
effect and composite pair superconductivity develop. Two 



vnqr = 



3eV zz Q 
2hI(2I - 1) 



Hz. 



(65) 



examples are shown - the black curve has T 2 > 7?, as ex- At the In(2) site, there is a nonzero asymmetric contri- 
bution, r\ — \V XX — V yy \/V zz , which can be independently 



pected for NpPdsAl 2 , so that the valence increases below T 2 , 
and then sharply decreases for T < T c . For the dashed gray 
curve, T* > T 2 * , as expected for CeMIns and a sharp in- 
crease in the valence is seen at T c . The sharp kinks at T* , T 2 
are artifacts of the mean field calculation, becoming smooth 
crossovers in real materials, while the superconducting kinks 
should be observable in real materials. 



moment / = 9/2, which results in a quadrupole moment, 
Q = 8.3 x l(T 29 m 2 , making them NQR activfP. The 
symmetry of the ground state doublet is T^I^S], whose 
angular dependence is given by, 



determined from experiment. Now that we have an ac- 
curate expression for the charge distribution of the two 
orbitals, we may calculate the electric field gradient at x 
associated with a charge q in orbital T at Rj : 



K r a(x,R, 



Aire 



OJx' 



■R,) 



3(a — a' 



l\2 



1 



rl\5 



r'|3 



(66) 



2^> cos 40 sin 2 9 sin 2£ (61) 
11 + 6 cos 2£ + 5 cos 20(1 + 2 cos 2£)] 



= — sin 
16 



where £ is a crystal-field parameter depending on the 
microscopic details that can be measured with inelastic 
neutron scattering, and we set £ s» .25 for CeCoIn 5 [2"5] . 
We take the symmetry of the second channel to be 1?2 = 
r 6 , whose angular dependence is, 



Summing over the eight neighboring Ce sites is sufficient 
to estimate the magnitude of the NQR shift, where we 
use the lattice constants of CeCoIn 5 , a = 4.6A, and c = 
7.4A[5]. The electric field gradients for a charge eqr in 
channel T are shown in Table I. 

For equal channel strengths, the total charge of the 
/-ion remains unity, and the increasing occupations of 
the empty and doubly occupied sites cause holes to build 
up with symmetry r^" and electrons with symmetry Tq, 
q = y6 = — 17+ ■ If we define q(T) to be the temper- 
ature dependent occupation of the empty/doubly occu- 
pied states, in the mean field q(T) — (6g) = (&|) will be 



proportional to T c 



T just below T c , and we can define, 
T r -T 



— [12 cos 20 + 5(3 + ce® 40)]. (62) 
32 



q(T) = q Q - 



(67) 



We can now use the real charge distributions of the two 
orbitals to estimate the magnitude of the electric field 
gradient, V a i, = —dE a /dxb, where 



where go is the ground state occupation of the 
empty/doubly occupied states. In terms of qo, the su- 
perconducting NQR frequency shift will be, 



pr(x)H$ r | 2 (0,0)i?(r) 2 . 
R(r) is the 4/ radial function, 



R(r) 



1 



fz 3 / 2 



96%/35 V 2a B 



exp 



Zr 
Aa B 



(63) 



(64) 



-7.6 g( / rc „ ^ kHz 



Aj/ 2 



NQR 



(T) 



-H?o 



T 

(T e -T) 



kHz 



(68) 



Even assuming a reasonably large 5% change in the 
single-occupancy (go = -05) for CeCoIns with T c = 2.3K 
will lead to a small shift in Vr} OR with a slope of sa 



as = -53A is the Bohr radius, and Z = 6vl0a.B /i"ce 
is adjusted so that the atomic radius is that of Ce 3+ . 
(r) = r Ce = 1.15 AgBJ. 



— .16, +.24 kHz/K, beginning precisely at T c . We could 
also consider the case of unequal channel strengths, where 
the Ce exchanges charge with the conduction electrons, 



12 



sitting in the In p-orbitals, however this term will be sim- 
ilarly small in any clean sample. If this shift could be dis- 
tinguished, it would be an unambiguous signal of compos- 
ite pairing. As this shift is quite small, Mossbauer spec- 
troscopy, which directly probes the /-ions may be a more 
likely technique to observe the development of a conden- 
sate quadrupole moment. Given that the crystal fields 
for NpPdsA^ are unknown, we can only make a rough 
estimate of the magnitude of the quadrupolar splitting, 
following [28] . The /-ion contribution to the splitting will 
originate from a redistribution of charge eq within the /- 
orbitals, and can be as large as 0.2q mm/s. The lattice 
contribution to the quadrupole splitting, which originates 
from the change in /-valence redistributing change within 
the lattice, will likely be a larger effect. 

VI. CONCLUSIONS 

Our two-channel Anderson model treatment has shown 
how composite pairing arises as the low energy conse- 
quence of valence fluctuations in two competing symme- 
try channels, which in turn manifests itself as a mixing 
of the empty and doubly occupied states, 

A sc oc(|0)(2|). (69) 

Composite pairing is primarily a local phenomena, where 
the pairing occurs within a single unit cell. The mixing 
is reminiscent of an intra-atomic antiferromagnetism, in- 
volving d-wave singlet formation between the \Ti) and 
|T2) states. It is the atomic physics of the /-ions, tuned 
by their local chemical environment that drives the super- 
conductivity. Such chemically driven d-wave pairing is a 
fascinating direction for exploring higher temperature su- 
perconductors in even more mixed valent 3c? materials, as 
the strength of the composite pairing increases monotoni- 
cally with increasing valence fluctuations, which accounts 
for the difference in transition temperatures between the 
cerium and the actinide 115 superconductors, and for the 
effects of pressure. 

The local nature of composite pairing should make 
it less sensitive to disorder than more conventional 
anisotropic pairing mechanisms. Such insensitivity is ob- 
served in the doped Ce-115 materials, where supercon- 
ductivity survives up to approximately 25% substitution 



on the Ce site^l, while neither Ce-site nor In-site dis- 
order behave according to Abrikosov-Gorkov, suggesting 
instead a percolative transition^. 

The redistribution of charge due to the mixing of 
empty and doubly occupied states provides a promising 
direction to experimentally test for composite pairing, 
which should appear as a sharp redistribution of charge 
associated with the superconducting transition. Both 
monopole (/-valence) and quadrupole (electric field gra- 
dients) charge effects should be observable, with core- 
level X-ray spectroscopy and the Mossbauer isomer shift 
or as a shift in the NQR frequency at surrounding nu- 
clei, respectively. We predict a shift with slope of or- 
der ±0.3kHz/K in the NQR frequency of In nuclei in 
CeCoIn 5 . 

Deriving these results in an exact, controlled mean field 
theory required the introduction of symplectic Hubbard 
operators, which maintain the time-reversal properties 
of SU(2) electrons in the large N limit. While our re- 
sults are obtained in the large- ./V limit, it should also be 
possible to use these Hubbard operators to develop a dy- 
namical mean field theory treatment of the two-channel 
Anderson lattice, enabling us to examine composite pair- 
ing for N = 2. 

In addition to the two-channel Anderson model, the de- 
velopment of symplectic-Hubbard operators allows a con- 
trolled treatment of the finite-?/ Anderson model, which 
is potentially useful as an impurity solver for dynamical 
mean field theory. We identify the finite-?/ model as a 
special case of our two channel model when the electron 
and hole fluctuations occur in the same symmetry chan- 
nel, Ti = r 2 . At first sight, this model appears to give s- 
wave superconductivity, however, the SU{2) constraint, 
Ai forbids any on-site pairing and completely kills the 
superconductivity, leaving only the simple Fermi liquid 
solution. We expect 1/N corrections to this mean field 
limit will differ from the SU (N) approach, and an inter- 
esting future direction is to use the Gaussian fluctuations 
to examine the charge fluctuation side peaks. 
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